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ABSTRACT 

This  paper  treats  the  motion  of  long,  unidirectional, 
continuous  fibers  through  a  Newtonian  resin.  There'  is  a 
consolidation  force  on  the  top  row  of  fibers  or  an  effective  force 
on  all  the  fibers  due  to  thei r  tension  in  a  f i i ament  wound 
cylinder.  An  asymptotic  analysis  for  high  fiber  volume  fraction  is 
presented  here.  The  results  of  the  present  rigorous  asymptotic 
solution  are  compared  to  the  results  of  a  previous  application  of 
the  1 ubr i c at i on  approx i mat i on . 
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INTRODUCTION 


I  . 


There  are  several  objectives  to  be  achieved  during  the 
consol  1  dat  i  on  of  a  composite  material.  The  elimination  of  voids, 
the  removal  of  air  and  excess  resin,  a  uniform  degree  of  cure,  a 
uniform  final  fiber  volume  and  conformation  to  a  specified  fiber 
orientation  are  all  necessary  in  the  final  consolidated  material. 
Ihe  consolidation  process  may  be  initiated  with  a  surface  or 
pressure  loading,  or  in  the  case  of  fiber  winding,  consolidation 
results  from  the  fiber  tension  due  to  the  winding.  Both  types  of 
loading  produce  resin  and  void  motion,  and  this  motion  depends  on 
the  resin  matrix  properties  (such  as  viscosity,  thermal 
conductivity,  cure  behavior,1,  the  initial  fiber  volume  fraction, 
fiber  orientation,  and  fiber  properties  (thermal  conductivity, 
diameter,  fiber  length). 

The?  word:  presented  here  is  concerned  with  the  resin  flow  and 
fiber  motion  in  a  fiber— resin  system  with  an  initially  high)  fiber 
volume  typical  of  composites  made  with  thermoplastic  resins.  We 
use  an  approximate  approach  to  the  analysis  of  the  resin  flow  which 
< 1 >  will  reveal  much  of  the  important  physical  phenomena  associated 
with  the  non— New ton i an  resins,  (2)  will  provide  reasonably  good 
quantitative  predictions,  and  (3)  will  indicate  key  simplifications 
to  the  full  boundary  value  problems  which  can  then  mate  a  mo)  e 
accurate  model  tractable. 

The  approximate  approach  r  ct^esent  r,  a.i  extension  of  a  model 


developed  by  Lindt  C1J  for  fiber-  motion  in  a  Newtonian  resin. 
Lindt  assumes  a  rectangular  array  of  fibers  arranged  in  parallel 
rows  and  constrained  to  vertical  motion  only.  Lindt  assumes  that, 
the  important  resin  flows  are  in  the  gaps  between  the  fiber  rows 
and  columns,  and  uses  lubrication  theory  to  treat  these  flows.  a 
correct  version  of  this  model  represents  the  leading  order  term  in 
an  asymptotic  solution  for  e  1  ,  where  e~  is  the  ratio  of  the  qap 
size  to  the  fiber  radius.  The  model  is  not  restricted  to  Newtonian 
fluids;  the  only  restriction  is  for  the  fiber  volume  to  be  large, 
corresponding  to  narrow  gaps  between  fibers.  As  thermoplastic 
prepreg  tape  for  filament  winding  has  little  or  no  excess  resin, 
this  asymptotic:  solution  will  be  reasonably  accurate  for 
t  h  er  mop  1  a  s  t  i  c:  resin  flows. 

We  begin  by  presenting  the  asymptotic  solution  for  a 
Newtonian  fluid  in  order  to  compare  our  results  with  the  results  of 
Lindt.  This  will  provide  a  check  on  this  work  and  on  Lindt  s  work, 
'he  solution  presented  here  is  a  rigorous  asymptotic  solution  tor 
t < < 1 ,  with  matching  between  the  subregions  in  each  narrow  gap.  in 

a  future  paper  ,  the  model  will  be  applied  to  non . Newtonian  fluids 

and  irregular  fiber  arrays.  In  section  2  the  flow  field  in  the 
narrow  gap  regions  will  be  determined,  with  the  forces  on  the 
fibers  due  to  these  flows  derived  in  section  3.  The  numerical 
analysis  and  consolidation  results  are  presented  in  section  4. 
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2.  PROBLEM  FORMULATION  FOR  NEWTONIAN  RESIN  FLOW 

We  consider  the  resin  and  -fiber  motion  in  a  system  o-f  n  rows 
o-f  moving  -fibers  suspended  in  a  high  viscositv  Newtonian  fluid. 
The  fiber  length  is  much  larger  than  the  fiber  radius,  and  the  flow 
along  the  fiber  length  is  negligible.  The  fibers  are  arranged  in 
columns,  with  the  separation  between  adjacent  columns  equal  to  260. 
The  adjacent  fiber  rows  are  separated  by  a  distance  2h1  (t.)  ,  whet  e  i 
denotes  the  row  of  the  fiber  above,  and  the  row  gap  is  a  function 
of  time.  The  fiber  columns  are  constrained  to  remain  in  columns 
with  the  same  gap  size  for  all  time  t.  The  geometry  of  the  fiber- 
resin  system  is  shown  in  figure  1.  All  fibers  are  subjected  to  an 
applied  force  due  to  the  fiber  tension  and  the  top  row  of  fibers  is 
subjected  to  a  distributed  surface  force  per  unit  length  due  to 
the  consolidation  process. 

The  Reynold's  number  tor  the  resin  is  Re=pU(_R/p,  where  f  and 
p  are  the  density  and  viscosity  of  the  resin,  while  Uf.  and  R  are 
the  character i sti c  velocity  and  fiber  radius.  For  consolidation 
with  polymer  ic  fluids,  Re  .  .1,  so  that  the  inertial  "force" 
associated  with  the  resin  acceleration  is  negligible  compared  to 
the  viscous  force  between  the  fiber  and  resin.  Since  the  density 
of  the  fibers  is  comparable  to  that  of  the  resin,  the  inertial 
"force"  associated  with  the  fiber  acceleration  is  also  negligible. 
Therefore,  there  is  cr  equi  1  i  briar,  of  -Fie  vertical  forces  for  c*o«2  ■ 


fiber  . 


£C 

vJ 


There  are  two  flow  regimes  for  the  given  geometry.  The  flow 
between  adjacent  rows,  caused  by  the  applied  loading  of  the  liber  s 
and  subsequent  vertical  motion  of  the  fiber  rows,  is  called  the 
squeezing  flow,  after  Lindt  LIT.  This  f 1 ow  is  the  source  of  fluid 
which  flows  through  the  gaps  between  the  adjacent  columns,  which  is 
called  the  normal  -flow,  again  after  Lindt  [13. 


h .  Squeeze  Flow  Between  Fiber  Rows 

For  the  squeeze  flow  between  fiber  rows,  the  top  fiber  (in 
row  i >  moves  downward  with  a  velocity  (t)  and  the  bottom  fiber 
(in  row  i+l)  moves  downward  with  a  velocity  V^+, (t>.  We  choose  a 
coordinate  system  with  an  origin  which  is  aligned  with  the  fiber 
centers  and  remains  midway  between  the  fiber  row  gap,  as  shown  in 
figure  2(a).  Thus  the  origin  moves  with  a  velocity  Lh  = 
--  ( V  +  j  )  /2  relative  to  a  fixed  coordinate  frame.  In  the  moving 
frame,  the  top  fiber  1  moves  downward  with  a  velocity 


while  the  bottom  fiber  l+i  moves  upward  with  a  velocity 


V.  -  V.  , 

1  l  f  1 


(lb ) 


The  squeezing  flow  between  the  i  and  i+1  fibers  is  symmetric 
about  the  >:=0  and  y=0  planes.  The  position  of  the  top  fiber 
surface  is  given  by 

R*-  =  (y— R— hj  (2> 


This  can  be  rewritten  as 


y  =  ti  1  +R— R 


<  l-x2/R2)  1 ' 2 


where  R  is  the  -fiber  radius  and  h^  is  halt  the  gap  size  between  the 
i  and  it  1  rows.  Equation  (3)  can  be  expanded  -for  :■:<  <R  to  qi  ve 


=  1  + 


2Rh 


The  equations  governing  the  inertiaiess  fluid  flow  are 


=  ~XP  + 


f  5  a--  c > 


where  v'L"= 


dx“  3y' 


and  u,p,p  are  the  velocity,  pressure  and  viscosity  of  the  fluid. 

We  consider  first  the  near  region  of  the  squeeze  region, 
where  is  comparable  to  h  and  both  and  h  x  are  much  less  than  R. 
For  !•;  comparable  to  hj,  the  fiber  surface  is  asymptotically  flat, 
and  the  surface  is  at  y/hz  =  l.  We  use  the  following  scalings  tor- 
region  1 


y,=  hiVl 


’  xl  “  'Vl 


^  /  a. .  tr 


* 

n  u 

i  y , 


U  =  ri  u 

X  1  X 


( 8a , b ) 


P1  = 


h  fJ  1 


,  t 


±  t* 


1 9  a  .  b  ) 


where  the  dimensionless  variables  are  denoted  by  an  aster  i  si 


The 


non-dimensional  equations  for  region  1, 


with 


the 
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asterisks  dropped,  are 
fp. 


( > 


o  = 


1 


+  V  U 


1 


1 


dx 

*Pl  2 

8y,  +  v  V 


V  •  a  =  0 
iv  >w 

2  3' 

where  V  =  - +  — 


(  1  Oa ) 

<  1  Ob ) 

(11) 


dx 


i>V  . 


1  ’1 

Due  to  the  symmetry  of  the  flow  and  geometry,  only  the  top  fiber 
row  i  need  be  considered  in  determining  the  squeeze  flow.  The 
boundary  conditions  are 


V  ’ 


u.  =  0  at  y  —  1 

*  1 


3u 


u  =  (.) 
V  , 


U  =  O  , 

::  1 


1 


9y 


J 

3u. 


0  at  y  -  0 


dx 


at  ^  -  0,  0  <y  ^  2 1 


( 1 2a  ,  b  > 


( 1 3a , b ) 


( 1 4a , b ) 


The  boundary  conditions  (1_-,14)  are  symmetry  conditions  on  the 
squeeze  flow  and  shear  stress. 

We  replace  each  dependent  variable  bv  a  laylor  series  with 
powers  of  ^  times  coef  f  l  cent  functions  of  Y\’  1  he  coeffi  cents  for 

the  different  powers  of  >;  j  decouple  and  are  governed  by  ordinary 
differential  equations  in  y^.  One  decoupled  problem  has  an 
inhomogeneous  boundary  condition  <l2a»  and  thus  has  a  non  _  et  o 
solution,  while  all  other  decoupled  problems  have  homogeneous 
boundary  conditions,  and  thus  have  only  zero  as  a  solution. 
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Therefore,  the  solution  represents  c  Taylor  series  in  x  ^  with  onl  y 
one  term  so  that  it  applies  at  all  x,.  The  only  non-zero  parts  of 
the  Taylor  series  are 


uw  =  x  U,  (y  > 

1 

u  =  V  ( y  4  ) 

y,  yi 


(lbs) 

(  1 5b  > 


R  ( y  ) 


(15c) 


When  the  solutions  (15)  are  substituted  into  (10,11)  and  the 
boundary  conditions  (12-14)  are  satisfied,  the  solutions  tor  region 
1  are 


(16a) 


(16b) 

(  16c ) 


where  8  is  a  constant  ambient  pressure  to  be  determined. 

When  :<>>hj  ,  the  curvature  of  the  surface  becomes  important. 

From  equation  (4)  ,  the  leading  order  curvature  term  is  0(1)  when  x 

1  /  2 

is  comparable  to  (Rhj  )  .  fhe  position  of  the  surface  is 

parabolic  in  x  for  region  2  where  hx  ; <  ;;=0  (Rhj  )1/ "  .  <R,  with 


h.  1  +  2Rh 

l  l 

The  appropriate  scalings  for  the  region  2  t i ow  variables  are 


(  1  7  > 


9 


=  <  Rh  ) 

o  1 


1  /2 


u 


* 

n  u 

1  W. 


M»'l 


P  = 


l _ 


eh 
l  l 


u  = 


"i  * 
-  Lt 

e  . 

l 


±  t* 


<  IB  a  ,b) 

i  1  9a  ,  b  ) 


c  2 9  a , b ) 


i  /  2 

where  E  i  “  i /R >  Thus  -for  region  2, 

compressed  by  the  -factor  e1  relative  to  y, 
scaling  -for  u„  -from  conservation  o-f  mass, 
dimensional  equations,  using  the  scalings 
(5,6),  are 


the  variable  is 

which  determines  the 
The  rescaled  non- 
i 1 8-20 )  i n  equdti ons 


0 


9P 


8u  du 

y,, 

<J  "  17T~  4  W7~ 


where  the  asterisks  have  been  dropped,  and  the  0 ( z  i  ) 
been  neglected.  The  boundary  conditions  are 


u  =  -1  ,  u  0  at  y„,=  t  ;  =  1  +  — 

y  . 


3uw 

ii  =  i)  ,  - —  —  0  at  y  k> 

y„:.  ‘  c»y ... 

u  -*  u  ,  u  a  u  as  ;  , , a  <_> 

y2  V1  *2  *1 


(  2  1  a  ) 


v  2  1  b  ) 


terms  have 


(  25a , b 1 


v  24a , h  > 
(  25a , b ) 
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Condition  (25)  is  a  matching  condition  between  region  2  and  region 

1 . 


Solving  (21,22)  with  the  boundary  conditions  (23-25)  yields 
the  solutions 


dp, 


P  .  =  — ~ 

2f  J 


( 26a ) 


(26b  ) 


( 26c  ) 


where  p represents  the  dimensionless  elevation  of  the  pressure  in 
the  squeeze  qap  above  the  ambient  pressure  for  >1.  The  pressure 
is  an  even  function  of  The  solutions  (26)  with  x -,=«.«  j  1  match 
the  inner  region  1  solutions  (16;  with  x  -xr,/  e  ,•  1,  and  this 
matching  determines  B=1  in  the  inner  region  1  pressure  (16c). 

the  fluid  squeezed  out  between  the  i  and  l+J  fibers  is  the 
source  of  fluid  flow  through  the  gaps  between  the  fiber  columns. 
It  is  this  normal  flow  which  is  considered  next. 

B.  Normal  Flaw  between  Fiber  Columns 

As  shown  in  figure  2(b) ,  the  adjacent  fiber  columns  are 
separated  by  a  distance  220,  where  6  remains  constant  tin  ouuhoul 
the  consolidation  process.  The  fibers  in  row  i  move  with  a 
velocity  vt  .  For  the  normal  flow,  we  use  a  coordinate  system  at 
rest  with  fiber  row  i,  with  the  origin  o+  the  coordinate  system 
midway  between  the  adjacent  columns.  The  surface  of  the  fiber  is 


given  by 


neg 1 ec  ted . 


Because  of  the  symmetry  o-f  the  flow,  only  one  side  of 


the  ith  fiber  column  need  be  considered.  Here,  we  treat  the  right 
side  of  the  column  cell.  The  boundary  conditions  on  the  normal 
flow  at  the  fiber  surface  are 

r~t 

u  =  0  ,  u  =  0  at  x=  f  (y)  =  1+  (33a,b) 

y  n  2  ’ 

n  n 

in  a  h-aine  of  reference  at  rest  with  the  fiber  row  l.  The  position 
of  the  surface  is  found  by  an  expansion  of  (27b)  using  the  scalings 
(28).  The  symmetry  condition  is  given  by 

du 

y 

-r -  ~  0  ,  u  =  0  at  ::=0  <34a,b) 

<7  y% 

n 

so  that  u  reaches  an  extremum  at  the  center  of  the  normal  qap. 

■  n 

There  is  an  additional  condition  on  the  normal  flow  which 
represents  a  constraint  on  the  amount  of  fluid  passing  through  the 
it.ii  column  qaps.  I  lie  fluid  which  is  squeezed  out  cd  the  region 
between  the  i  and  i+l  fiber  rows  is  forced  into  the  l  column  gaps, 
and  hence  flows  upward.  In  the  region  between  the  i+l  and  i +2 
rows,  fluid  is  also  lost  due  to  squeezing  and  is  forced  into  the 
l  a  1  column  gaps  arid  upwards  through  the  i  normal  gaps,  the  i-l 
normal  gaps,  and  so  on  until  the  fluid  reaches  the  top  fiber  row. 
The  total  volume  occupied  by  the  resin-fiber  columns  changes  as  the 
fluid  is  removed  frum  the  squeeze  region,  while  the  size  of  the 
column  gaps  remain  constant.  The  change  in  volume  of  the  column 
front  row  l  to  row  n  (the  static  bottom  fiber  row)  in  a  time  At  is 
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equal  to  the  volume  o-f  fluid  which  is  squeezed  and  then  passes 
through  the  column  gaps.  Thus  ,  the  volume  flow  rate  per  unit 
length  passing  through  the  ith  column  gap,  is  given  by 


Q  = 
1 


L 

o 


V.  At 

l 

At 


where  2<R+&0)  is  the  horizontal  size  of  the  column  cell  about  fiber 
l  and  V j At  is  the  total  change  in  the  vertical  size  of  the  l  to  n 
column  cells.  Then  G'r  represents  the  accumulated  flow  due  to  the 
squeezing  between  the  fiber  rows  from  i  to  n.  The  above  expression 
for  GL  can  be  simplified  and  rewritten  as 

Q  =  2  I R  +  L  I V  ( 35 ) 

i  !  o  |  i 


It  should  be  noted  that  the  total  volume  flow  rate  per  unit  length 
can  be  calculated  from  the  top  fiber  row  velocity  with  i-l;  this 
determines  the  net  loss  of  resin  in  the  system. 

The  volume  flow  rate  per  unit  length  Q1  corresponds  to  the 
flux  of  fluid  through  the  i th  column  gap.  Thus,  using  the  scalings 
<..'■8-30)  ,  the  conservation  of  mass  condition  is 


f  ( y  > 
n 

u  d ;  : 

vn 

-  o 


Q  £ 

i  n 


V  C 
1  o 


1  36 ) 


F  or  s i mp 1 icit/ 


we  use  the  notation 
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Q.  £ 

1  n 

qi"  2V. 6 

l  o 


£  + 
n 


1 


£ 

n 


<  3  7 ) 


Equations  (31,32)  are  solved,  using  the  boundary  conditions 
(33,34)  and  conservation  of  mass  (36)  ,  -for  the  normal  gap  -flow. 
The  solutions  are 


u  = 
"  n 


.  d*"p 
1  *n 


dy' 


.d  l 

3  I 


dP, 


+  yf 


n  dy 


u  = 
yn 


1  dpn 


2  dy 


dPr 


-q. 


dy 


<  38a ) 


(38b) 


( 38c ) 


r-rorn  equation  (38c)  it  is  seen  that  there  is  a  pressure  drop  across 
the  normal  gap  which  drives  the  flow  upward  toward  the  surface  of 
the  system,  and  that  this  pressure  depends  on  the  volume  flow  rate. 
As  yaO,  the  solution  (38)  becomes  simple  Poiseuille  flow  between 
parallel  walls  at  • .  -  .1.  J  ,  which  is  the  correct  solution  for  the 
normal  flow  when  the  dimensional  ;;  and  y  are  both  comparable  to  b0. 
While  the  squeezing  flow  requires  a  special  solution  in  the 
interior  region  1,  the  normal  flow  does  not  require  an  inner 
region,  as  previously  stated. 


FORCE  BALANCE  ON  THE  FIBERS 


Since  the  inertial  "force"  due  to  the  fiber  accelerations  is 
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negligible  compared  to  the  viscous  forces  of  the  resin  on  the 
fiber,  the  sum  of  the  forces  on  the  fibers  must  be  zero.  the 
forces  on  the  fiber  have  two  sources: (1)  the  externally  applied 
force  due  to  a  fiber  tension  or  distributed  surface  force  per  unit 
length,  and  (2)  the  force  the  fluid  exerts  on  the  fibers.  The 
externally  applied  force  per  unit  length  is  denoted  as 


F1 
/v>  ex 


t 


(39) 


where  the  applied  load  acts  in  the  -y  direction  with  a  magnitude  of 
I  .  at.  each  fiber  row. 

The  forces  the  fluid  exerts  on  the  i th  fiber  surface  are 
found  using  the  relation 


dF1  =  a  *  n .  d  A  (40) 
J  J  k  k 

where  a1  is  the  stress  tensor  for  the  fluid  due  to  flow  about  the 
ith  fiber.  Equation  (40)  is  the  incremental  force  in  the  ith 
direction  which  acts  on  an  incremental  surface  area  element  dA 
which  has  a  unit  vector  rq  , . 

A.  Force  due  to  the  squeezing  flows 

There  are  two  parts  to  the  force  on  the  l th  fiber  row  due 
to  the  squeezing  flows.  The  first  part,  is  the  force  on  the  ith 
fiber  due  to  the  squeeze  flow  between  the  i  and  i+1  fiber  rows 
acting  on  the  bottom  of  the  ith  fiber.  The  second  part  of  the 
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•force  is  due  to  the  squeeze  flow  between  the  i-1  and  i  fiber  rows, 
acting  at  the  top  of  the  ith  fiber. 

The  first  step  to  finding  the  forces  due  to  the  squeezing 
flow  is  to  find  the  normal  vector.  The  normal  vector  to  a  surface 
is  given  by 


A 

n  = 

where  G=0  is  the  equation  of  the  surface.  Using  equation  (3) , 
surface  equation  for  the  squeeze  flow  cell  is  given  by 


(41 ) 


the 


G  =  R 


h 

l 


0 


(42) 


The  normal  vector  to  the  top  surface  in  the  ith  squeeze  region  is 


*  A  ( R  —  )  a 

n  ~  R  X  ~  - R -  y 

A  A 

where  x  and  y  ar  e  unit  vectors.  by  substituting  the  scalings, 
into  the  expression  (43a)  ,  the  normal  vector  becomes 

A  *A  A 

n  =  sc  x  x  -  y 
12 


(43a) 

(18) 


(43b) 


where  the  0  (  e  1  )  term?;-  have  been  neglected  in  order  to 
consistent  with  previous  approximations.  The  area  element 
the  surface  is  found  by  a  similar  analysis.  Tor  a  circular 


r emai n 
dA  for 
surf  ace 


dA=Rd©  ,  where  ©=tan 


(y,_— R— h  ) 

2  1 


ae  .  ae  . 

aZ,d* 2  *  sy,dy2 


<  44  > 


and  d 6= 
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Ignoring  0 ( e x ) ^  terms,  the  area  element  (44)  is  given  by 

d  A=d;<  2  (45) 

Thus,  with  (40, 43b, 45)  ,  the  -forces  on  the  ith  -fiber  row  due 
to  the  squeezing  flow  are 


( 46a ) 


( 46b  ) 


where  the  limits  are  the  asymptotic  limits  of  the  surface,  and  the 
asterisks  denote  the  non-dimensional  variables. 

For  an  incompressible,  Newtonian  fluid  the  stress  is  given  by 


a 


jk 


( 47 ) 


The  scalings  (18,19,20a)  and  the  squeeze  flow  solutions  (26) 
are  substituted  into  (47).  The  integrand  in  (46a)  is  an  odd 
function  of  and  thus  gives  a  zero  net  contribution  to  the  x  force 
(46a).  This  is  seen  by  examining  the  solutions  (26):  the  pressure 
is  an  even  function  in  as  is  the?  normal  stress  (  >)u  ,  /  dx  0 )  ,  while 
the  shear  stress  {  du  /dy-,)  is  odd  in 

There  is  a  force  in  the  y-di r ect l on  due  to  the  pressure. 
With  the  pressure  given  by  equation  ( 26c ) ,  equation  (46b)  becomes 
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F 


i 

v 

sq ,  bot 


3rrnr|. 


<  46a ) 


A 

This  gives  the  -force  on  the  bottom  o-f  the  ith  -fiber  row  in  the  y- 
direction  due  to  the  pressure  generated  by  the  squeezing  flow 
between  the  i  and  i  +  l  fiber  rows.  The  force  on  the  top  of  the  ith 
fiber  row  is  directed  downward,  and  is  given  by 


3ttj.o-|.  . 

F  =  -  - i— -  (48b) 

Ysq  ,  top  2 /2  e  1 

This  force  is  due  to  the  squeezing  flow  between  the  l-l  and  i  fiber 
rows.  We  have  used  the  pressure  (26c)  throughout  the  squeeze  gap 
region,  but  there  is  a  very  small  region  with  >;=0(h^)  where  the 
pressure  is  given  by  equation  (16c)  with  B=l.  The  pressure  in  the 
region  x=0(h1)  is  slightly  less  lhan  the  pressure  predicted  by 
(26c).  However,  the  deviation  is  0  v  e  L and  extends  over  a  small, 
□  <e1)  portion  of  the  surface,  so  that  the  error  is  0(£j)'~'  smaller 
than  the  forces  in  equations  (48). 

B.  Force  due  to  the  normal  flow 

To  calculate  the  force  due  to  the  normal  flow,  we  proceed  as 
we  did  for  the  squeezing  flow.  Usinq  the  expression  (27b),  the 
equation  for  the  right,  fiber  surface  in  the  normal  gap  is 
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r>  o  i  / 

G  —  L  +R—  (R'4’—  y  - )  O  (49) 

O  7 

The  unit  normal  vector,  -found  using  the  equations  (28,41,49),  xs 
given  by 

A  A  *  A 

n  =  -x  +c  y  y  (50) 

neglecting  □ ( e  ) ^  terms,  and  the  asterisks  denote  the  non- 
dimensional  variables.  The  area  element  is  dA=dy.  Using  equations 
(40,50),  the  forces  due  to  the  normal  flow  at  the  fiber  surface  are 


r  '  •• 

-aww 

* 

+  e  y  a 

n  y  I 

(51a) 


F  = 

Y  n 
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*  1 
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£  V 

a  i 

l  yx 

n ' 

yy  1 

d  y 


(51b) 


where  the 

1  i  m  i  1 5 

are 

t  h  e  a  s  y  m  p  t  o  1 1  c. 

limits  for 

the 

r eg i on . 

The 

seal ings  ( 

28-3o a ) 

and 

the  solutions 

(38)  for  the 

normal 

f  1  ow 

are 

subst i tuted 

into 

equa 

t l ons  (47,51) . 

The  result 

f  or 

the 

force 

l  n 

the  x-di recti  on  is 
mV  f  + 

F1  =  p  dy  (52) 

w  |  n 

n  e 


where  the  integrand  is  in  terms  o+  the  non-dimensional  variables. 


and  the  asterisks  have  now  been  dropped.  By  examining  the  solution 


(38c) ,  the  integrand  of  (51a>  is  shown  to  be  an  odd  function  of  y. 


Thus,  the  force  P\,  due  to  the  normal  flow  is  zero. 

''n 

A 

The  force  in  the  y-di recti  on  due  to  the  normal  flow  is 


"VP,-  ~  37 

n  o  > 


where  the  integrand  of  (53)  is  in  terms  of  the  non-dimensional 
variables  (28-30),  and  0  (  c  )  ^  terms  have  been  neglected.  This 
force  (53)  is  produced  by  the  pressure  distribution  and  shear 
stress  at  the  fiber  surface. 

Before  evaluating  equation  (53),  the  pressure  distribution  in 
the  normal  gap  is  determined  by  integrating  the  equation  (38c)  to 
obtain 

— -  +  —  tan-1  — £■  l  (54) 

4  I  2  +  y"-|  4  /2  .  /2  J  J 

The  first  two  terms  in  the  expression  (54)  approach  zero  as  ya±'», 
but  the  last  term  approaches  two  different  constant  pressures  as 
y-*±*>,  so  that 

9  7Tq  i 

Ap  =  p  ( y  -»•>:•)  -  p  ( y  -*•—•»)  -  -  - — : —  (55) 
n  n  4/d 

This  Ap  is  the  pressure  drop  needed  to  drive  the  flow  through 
the  normal  gap  between  the  adjacent  fibers  in  the  ith  row. 
Therefore,  the  ambient  pressure  above  the  l th  fiber  is  less  than 
the  ambient  pressure  below  the  ith  fiber.  The  solutions  (26,38,54) 
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give  the  velocity  and  pressure  in  the  squeeze  and  normal  gaps,  gaps 
which  represent  a  small,  0(e)  -fraction  of  the  total  fiber  surface. 
The  pressure  drop  (55)  gives  a  pressure  over  the  entire  bottom  of 
the  fiber,  a  pressure  which  is  greater  than  the  pressure  over  the 
entire  top  of  the  fiber  by  a  value  j  Apn  |  .  This  pressure 
difference  produces  a  net  upward  force  on  the  fiber.  We  cannot 
obtain  this  force  by  introducing  the  expression  (54)  into  equation 
(53) ,  as  the  unit  normal  vector  (50)  is  valid  only  in  the  normal 
qap  region  where  the  fiber  locally  appears  to  be  a  parabola.  If  we 
were  to  introduce  the  pressure  (54)  into  the  force  expression  (53), 
we  would  obtain  an  infinite  force,  because  there  is  a  pressure 
difference  across  a  parabolic  surface,  which  has  an  infinite  area. 

To  find  the  total  vertical  force  due  to  the  normal  flow,  we 
first  consider-  the  vertical  force  due  to  the  change  in  the 
d i men si  anal  ambient  pressur  e 


on  the  top  and  the  bottom  halves,  respectively,  ot  the  fiber  with 
the  true  circular  cross  section.  This  gives  a  pressure  force  per 
unit  fiber  length. 


9rTq  pV 

i  l 


(57) 


Since  the  force  due  to  the  pressures  (56)  has  been  accounted  for  in 
the  pressure  force  (57),  the  expression  (53)  is  replaced  by  the 
expressi on 
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(58) 


for  the  force  per  unit  length  on  the  left  side  of  the  ith  fiber, 
due  to  the  local  parts  of  the  normal  gap  flow.  The  integral  (58) 
with  the  expressions  (38,54,55)  is  well  defined  and  gives 
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( 59 ) 


By  symmetry,  the  force  per  unit  length  on  the  right  side  of  the  ith 
fiber  is  also  given  by  the  expression  (59).  The  total  normal  force 
due  to  the  ambient  pressure  and  local  parts  of  the  normal  flow  is 
given  by 

9  jt  pV 

F1  -  F 1  4  21- 1  =  - i*- 

y  y  y 

n,tot  p  n  ~zyy.z^ 


n 


<  60 ) 


As  the  inertial  force  associated  with  the  fiber  accel erat l ons 
is  negligible,  the  downward  force  per  unit  length  (39)  must  equal 
the  upward  forces  due  to  the  squeeze  and  normal  flows,  where 


=  F 


.  l 


+  F 


+  F 


sq , bot 


sq , top 


n  ,  tot 


(61) 


Introducing  equations  (1,37,48,60),  the  force  balance  (61)  becomes 
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Equation  (62)  is  the  force  balance  for  a  system  of  n— moving 
fibers  with  squeeze  gaps  2h  which  satisfy  the  condition 

I/O  *  /  ■  . 

h  i  •  ( Rhx  )  A ^  <  <E  ,  and  normal  gaps  26Q  which  satisfy  L  ■  (R&  )  1 ^  •  R. 

The  force  balance  given  in  (62)  is  applicable  for  all 
interior  fibers,  i =2  to  n-1.  The  nth  fiber  is  considered  to  be 
stationary;  the  bottomi  fiber  is  constrained  by  a  solid  surface  and 
so  remains  motionless,  with  the  fibers  above  it  moving  toward  the 
nth  fiber.  The?  top  fiber  i~l  has  no  squeezing  flow  due  to  a  fiber 
above  it.  Thus,  for  i  =  1  ,  the  force  balance  • o2 )  is  modified  to 
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4.  NUMERICAL  CALCULATIONS.  AND  CONSUL  1  DAT  1  ON  RESULTS 


The  consolidation  behavior  of  the  set  of  n-movinq  fibers  is 
determined  by  solving  the  set  of  n-1  simultaneous  equations 
(62,63) ,  with  the  condition  that  the  i =n  fiber  is  stationary,  tor 
the  fiber  velocities,  fiber  displacements  and  volume  flow  rate  per- 
unit  length.  To  facilitate  the  calculation,  the  following  scalings 
are  used  in  (62,63) , 

H  r*. 

T=  7  T  ,  t  =  (o4a.br 

1  C  1  cl 
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j  1  /  2 

( 644 , g ) 

where  the  asterisks  denote  the  non-dimensional  quantites. 

The  value  of  T  is  taken  to  be  the  magnitude  of  the  tension 
in  the  -fibers  or  of  the  distributed  surface  force  per  unit  length. 
The  char acter 1 st 1 c  time  is  determined  using  the  fluid  viscosity, 
fiber  radius  and  characteristic  applied  force  T  .  Since  the 
characteristic  time  depends  on  the  fluid  viscosity,  a  change  in  the 
viscosit .  will  change  the  value  of  t  ,  and  the  char acter i st i c 
velocity,  R/t  .  Similarly,  a  change  in  the  value  of  the 
char  a  ter  1  st  1  c.  applied  force  will  change  t  as  well  as  the? 
charac ter istic  vel oc 1 1  y. 

With  the  scalings  i64>,  the  force  balance  equations  162,63) 
arp  rewritten  as 
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As  fluid  is  squeezed  out  from  between  the  fiber  rows,  the  oap 


between  the  rows  will  change.  Vo  take  this  into 
squeeze  gap  2h  is  updated  at  each  time  step,  where 

2h  (t+At>  =  2h  it)  +At(V  4  -  V  ) 

1  1  1  +  1  1 

or  using  the  scalings  (64; 

*  *  *  *-*  *  ■*  #■ 

h  ( t  +  At  )  =  h  ( t  )  +  At  (V  -  V  )  /  2 

l  l  l  + 1  l 

The  equations  <65,661,  along  with  tne  zero  velocity  condition 

on  the  nth  -fiber,  are  -formed  into  a  set  of  simultaneous  equations 

of  the  form  L  v*  =  T*.  The  fiber  velocities  are  found  usinq  a 

1 1  i  J  - 

routine  tor  solving  sets  of  simultaneous  linear  algebraic  equations 
[21.  Once  the  velocities  are  known,  the  squeeze  gaps  2h*,  using 
equation  (67),  are  updated  and  the  next  time  step  is  performed 
until  t  *- 1 

(Tic*. 

A.  Consolidation  from  the  top:  Lindt's  problem 

We  con  rider  first  the-  calculation  corresponding  to  one  of 
Lindt  s  systems  for  n=12  -libers.  I  he  top  fiber  row  <i--J)  is  loaded 
by  a  distributed  surface  force  per  unit  length  ,  I  ^ -O .  <X>5N/ m ,  and 
there  are-  no  other  applied  forces  on  the  interior  fibers  i=2-12. 
Tf  fluid  viscosity  is  Julg/rn— s  arid  the  fiber  radius,  R,  is  taken 

_ r. 

to  be  10  “'m.  ffie  initial  gaps  are  given  by 

h.  <  t  *—0 )  -  L  ~  R  ~lu  m  (6b) 

1  o 

The  diaracter istic  time,  t  ,  i s  0.02  seconds.  Half  the 

dimensionless  distance  between  the  fiber  rows,  h* ,  for  (i=l— 5)  and 


account ,  the 


(67) 
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the  dimensionless  fiber  velocities,  V*,  for  <i=l-5>,  are  shown  in 
figure  3(a)  and  3(b)  respectively.  The  volume  flow  rate  per  unit 
length  may  be  calculated  directly  from  the  fiber  velocities  using 
(35)  and  the  scalings  (64).  The  dimensional  flow  rate  per  unit 
length  is  given  by 

Q.  =  10"e(l  +  &*)V*  nws  (69) 

1  o  1 


The  dimensionless  normal  force  (60)  and  dimensionless  squ^e:e  -force 
(48)  for  the  i=l  fiber  is  shown  in  figures  4(a),  and  figure  4(b) 
shows  the  dimensionless  normal  and  squeeze  forces  for  the  i=2-4 
fibers. 


The  results  for  this  problem  can  be  intrepreted  as  a  type  of 
consolidation  wave'.  Consider  the  i=l  fiber  for  t£0,  when  the 
velocity  of  the  i ~2  fiber  is  close  to  zero.  From  (65),  the  ratio 
of  the  normal  force  to  the  squeeze  force, 


(  70) 


For  en=l,  initially  e^l,  the  normal  force  is  24  times  larger  than 
the  squeeze  force.  The  l = 1  fiber  is  then  in  a  free  fall  due  to  the 
applied  force,  with  the  resistive  force  due  to  the  normal  flow 
only;  the  interaction  between  fiber  1  and  fiber  2  due  to  the 
squeezing  flow  is  very  small.  The  velocity  V*  is  then  controlled 
by  the  normal  flow,  and  the  velocity  can  be  determined  by  the 
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balance  between  the  applied  force  and  the  normal  force,  with  1*=1 
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then  V*=0.025  for  £^=1.  The  squeeze  force  acting  on  the  bottom  of 

the  i=l  fiber  and  top  of  the  i =2  fiber  is  small,  but  it  does 

influence  the  i =2  fiber.  As  the  i=l  fiber  moves  downward,  the 

squeeze  force  increases  in  maqnitude,  since  F  Thus,  the 

sq  l 

force  balance  requires  that  the  normal  force  acting  on  the  i=l 
fiber  decrease  as  the  squeeze  force  increases.  From  equation  (60) 
it  is  seen  that  Fn<>:VF,  so  the  velocity  of  the  i  =  l  fiber  must 
decrease  to  decrease  the  normal  force.  The  squeeze  force  acting  on 
the  i=l  and  i=2  fibers  begins  to  increase  the  velocity  of  the  l =2 
fiber  from  zero,  in  turn  increasing  the  normal  force  on  the  i =2 
fiber,  while  the  i =2  fiber  velocity  is  less  than  the  1=1  fiber 
velocity.  When  the  squeeze  qap  between  the  i=l  and  l =2  fibers 
becomes  smal 1 ,  the  squeeze  force  acting  on  the  1=1  fiber  begins  to 
dominate.  This  occurs  at  t  >c2 . 4  seconds  for'  the  i  =  l,2  fibers.  At 
thus  time,  the  squeeze  force  acting  on  the  l  =2  fiber  is  mainl,  due 
to  the  interaction  between  the  i=2  and  i=l  fibers,  with  the 
interaction  between  the  l  =2  and  i  =3  fibers  very  small,  and  the  l  =2. 
fiber  almost  stationary. 

The  applied  load  T*  is  distributed  between  the  i=l  and  l =2 
fibers,  after  t^2.4  seconds.  The  normal  force  on  the  two  fibers  is 
equul  ,  but  only  halt  the  value  o-t  the  initial  normal  force  on  the 
i=l  fiber  only.  The  squeeze  force  acting  on  the  two  fibers  is 
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equal  in  magnitude  but  opposite  in  direction.  The  two  fibers  are 
then  carrying  the  same  load,  and  begin  to  move  downward  with  the 
same  velocity.  This  velocity  is  roughly  one-half  the  free-fall 
velocity  for  one  fiber,  as  now  the  applied  load  has  been  split  into 
two.  These  two  fibers  have  formed  a  consolidation  front. 

As  the  top  two  fibers  move  downward,  the  i=3  fiber  begins  to 
be  influenced  by  the  decrease  in  the  squeeze  gap,  2h*,  between  the 
i=2  and  i=3  fibers.  The  squeeze  forces  acting  on  the  bottom  part 
of  the  i=2  fiber  and  the  top  part  of  the  i =3  fiber  increases  in 
magnitude  as  the  squeeze  gap  decreases.  As  yet,  the  i =4  fiber  is 
stationary  and  does  not  contribute  to  the  overall  force  balance. 
Since  the  squeeze  force  on  the  i~3  fiber  is  non-zero,  the  fiber 
begins  to  move,  slowly  at  first,  but  increasing  as  the  i =2  fiber 
approaches.  The  velocity  of  the  i=3  fiber  produces  a  normal  flow, 
and  hence  a  normal  force,  which  opposes  the  squeeze  force  acting  on 
the  top  of  the  i=3  fiber  (with  no  squeeze  force  acting  on  the 
bottom  of  the  3=3  fiber)  and  gives  a  balance  of  the  forces.  The 
total  squeeze  force  acting  on  the  l =2  fiber  decreases  with  the 
addition  of  the  opposing  squeeze  force  on  the  bottom  of  the  i=2 
fiber.  This  decreases  the  velocity  of  the  i =2  fiber  further,  and 
through  the  force  interaction,  the  i=l  velocity'  is  also  reduced. 
Because  the  l =3  fiber  has  no  interaction  with  the  1=4  fiber  as  yet, 
the  i =3  fiber  velocity  steadily  increases,  although  remaining  less 
than  the  velocities  of  the  1=1,2  fibers. 

When  the  squeeze  gap  between  the  i =2  and  l =3  fibers  becomes 
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small,  the  total  squeeze  -force  acting  on  the  i  =3  and  i  =2  -fibers  are 
equal.  This  occurs  at  t%7  seconds.  Further,  -force  balance 
requires  that  the  normal  -forces  are  also  equal.  The  normal  force 
acting  on  the  i=l-3  fibers  is  one-third  of  the  initial  normal  force 
acting  on  the  i=l  fiber.  The  applied  force  at  t%7  seconds  has  then 
been  distributed  among  the  three  fibers,  forming  a  three  fiber 
consol i dat i on  front,  which  moves  with  a  velocity  which  is  one- third 
the  free-fall  velocity  for  one  fiber. 

This  consolidation  process  continues  as  the  consolidation 
front  moves  downward  toward  the  fiber  directly  below  the  front. 
Unce  the  squeeze  gap  becomes  sufficently  small  tor  the  total 
squeeze  force  on  the  bottom  fiber  in  the  front  to  equal  the  squeeze 
force  acting  on  the  top  of  the  fiber  directly  below  the  front,  the 
fiber  below  becomes  a  part  of  the  consolidation  front.  This 
process  simply  propagates  downward  row  by  row,  with  a  velocity 
which  decreases  as  each  row  is  added.  The  front  starts  to  move 
with  a  free-fall  velocity  determined  by  the  balance  between  the 
normal  force  and  an  applied  load  distributed  among  the  G-fibers  in 
the  front.  This  front  free-fall  velocity  is 


where  G  is  the  number  of  fibers  in  the  consolidation  front. 

The  ley  to  the  consol l dat i on  is  in  the  normal  flow.  From 
(60)  it  is  seen  that  the  normal  force  is  strongly  dependent  on  the 
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size  of  the  normal  gap.  For  large  normal  gaps,  the  fluid  from  the 
squeezing  flow  encounters  little  resistance  as  xt  flows  into  the 
normal  gap  region,  and  so  the  consolidation  front  can  propz.gate 

rapidly.  However,  as  the  normal  gap  is  decreased  in  size,  the 

resistance  term,  here  defined  as 

e4  +  l  r 

(73) 

£ 

n 

becomes  large,  and  so  the  free  fall  velocity  of  the  front  (72)  is 
reduced.  This  indicates  that  the  fluid  in  the  squeeze  region 
cannot  penetrate  into  the  normal  gap,  and  so  consolidation  is 
inhibited.  Figures  5(a-c)  show  the  behavior  of  a  system  with  the 

normal  gap  reduced  by  half,  where  6  =0.5R.  This  case  has  the  same 

loading,  viscosity,  fiber  radius  and  character l st l c  time  and 
initial  squeeze  gap  (h(t— 0>=R)  as  the  previous  case,  and  the  fiber 
behavior  conforms  more  closely  to  the  tree  fall  analysis  already 
discussed.  This  is  due  to  the  increase  in  the  initial  normal  force 
caused  by  the  decrease  in  the  normal  ciap ,  where  the  ratio  of  the 


normal 

force  to  the 

squeeze  forct 

3  (70)  is  increased 

by 

about  a 

f  actor 

of  3.2.  The 

i  =  1 

f i ber  i s 

in  free  fall,  with 

a 

vel oc i ty 

rough  1 y 

one-third  of 

the 

case  for 

£/o=R;  this  decrease 

1  s 

directl y 

related  to  the  increase  in  the  normal  gap  resistance  to  the  flow 
from  the  squeezing  action  between  fiber  rows  1  and  2.  Since  the 
velocities  are  reduced,  the  time  required  for  the  i=l  fiber  to 
close  the  squeeze  gap  suff icently  for  the  squeeze  force  to  dominate 
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increases.  Once  the  i  =  l,2  -fibers  -form  the  -front  <  at  t£7  seconds; 
the  -front's  velocity  is  half  of  the  initial  free  fall  velocity,  as 
determined  by  (72)  with  C=2.  Again,  the  fibers  are  moving  slowly 
in  a  free  fall,  until  the  next  squeeze  gap  is  small.  The  value  of 
the  size  of  the  squeeze  gap  at  which  the  normal  and  squeeze  forces 
are  equal,  is  also  reduced  as  the  normal  gap  size  is  reduced,  as 
indicated  in  (70). 

If  the  normal  gap  is  reduced  more,  the  behavior  of  the 
consolidation  would  be  similar  to  the  first  two  cases,  but  with  an 
increase  in  the  time  for  the  fibers  to  consolidate,  and  a  stronger 
conformation  to  the  free-fall  behavior  previously  described. 

B.  Consolidation  with  uniform  fiber  tension 

The  next  case  considered  is  one  in  which  the  n-1  fibers  have 
a  uniform  tension  applied.  The  fiber  radius  has  a  value  of 

_ rzr 

R=10  with  h  <  t=0)  =  £>  -R ,  fluid  viscosity  is  lOkg/m-s  and  the 
applied  force  due  to  the  tension  of  winding  on  a  curve  is  0.005N/m, 
such  that  1 1  for  i=l  to  n-1.  Hal 1  the  dimensionless  distance 
between  fiber  rows,  h*  ,  and  the;  dimensionless  fiber  velocities,  V* , 
are  shown  in  figures  6(a)  and  6(b)  respectively,  with  the 
dimensionless  normal  and  squeeze  forces  shown  in  figure  6(c). 

The  consolidation  under  uniform  fiber  tension  is  similar  to 
the  consolidation  for  uniform  top  loading.  However,  with  a  uniform 
tension  on  all  of  the  fiber  rows,  the  consolidation  starts  at  the 


i = 1 1  fiber  and  pr ogresses  upward.  initially,  ail  n-1  fibers  are  i n 


a  uni-form  -free-fall  (  i=n  -fiber  is  static),  with  a  -free-fall 
velocity  determined  by  (71)  with  f*  replacing  T*.  It  is  the  i=ll 
fiber  which  experiences  a  squeezing  for^e  due  to  the  presence  of 
the  static  i=12  fiber.  The  squeeze  force  acts  on  the  bottom  of  the 
i=ll  fiber,  with  the  squeezing  force  between  the  i=10  and  i=il 
fibers  initially  zero  and  remaining  small.  As  shown  in  figure 
6(c),  the  behavior  of  the  normal  and  squeeze  forces  for  the  i = 1 1 
fiber  is  analagous  to  the  force  behavior  of  the  i=l  fiber  in  the 
uniform  top  loading  case.  As  the  i = 1 1  fiber  moves  downward  toward 
the  static  i-12  fiber,  the  squeeze  force  increases  as  the  squeeze 
gap  decreases.  This  necessitates  the  decrease  of  the  i=ll  fiber 
velocity,  to  decrease  the  normal  force  and  maintain  the  force 
balance.  However,  unlike  the  top  loading  case,  as  the  i = 1 1  fiber 
approaches  the  1  =  12  fiber,  the  squeeze  force  acting  on  the  l = 1 2 
fiber  does  not  produce  the  associated  motion  of  the  i=12  fiber,  due 
to  the  boundary  condition  that  the  i - 1 2  or  i =n  fiber  remain  static. 
Thus,  the  gap  between  the  i=ll  and  i=12  fibers  continues  to 
decrease,  with  the  only  resistance  to  the  motion  produced  by  the 
increasing  upward  directed  squeeze  force  on  the  bottom  of  the?  i=ll 
fiber.  The  ratio  of  the  normal  force  to  the  squeeze  force  for  the 
i=ll  fiber  under  uniform  tension  is  the  same  as  that  given  in 
equation  (70),  with  i=l  replaced  by  i=ll.  When  the  squeeze  gap, 
2h*j,  is  .24  (  at  t%1.6  seconds) ,  the  normal  and  squeeze  forces 
acting  on  the  i = 1 1  fiber  are  equal.  As  the  squeeze  gap  continues 
to  decrease,  the  normal  force  decreases,  while  the  squeeze  force 


increases  and  the  squeeze  gap  asymptotically  approaches  zero 
extent,  with  a  velocity  which  also  tends  to  zero.  The  normal  -force 
approaches  zero,  and  the  large  pressure  in  the  squeeze  gap  between 
the  i  =  ll  and  i  =  12  fibers  produces  a  sufficent  force  to  balance  the 
uniform  tension.  The  i  =  ll  and  i  =  12  fibers  form  the  consol  1 dat 1  on 
front . 

The  other  fibers  are  completely  unaffected  by  the  interaction 
between  the  i  =  ll  and  i=12  fibers,  and  the  i=l-10  fibers  continue  to 
move  in  a  free-fall  with  the  velocity  given  by  (71).  However,  once 
the  i=ll  fiber  has  virtually  stopped,  the  i=10  fiber  closes  the 
squeeze  gap  with  the  now  static  i  —  1 1  and  i  =  12  fibers.  The  decrease 
in  the  squeeze  gap  between  the  i=10  and  i=ll  fiber  increases  the 
squeeze  force  from  zero,  and  this  force  acts  on  the  bottom  of  the 
i=10  fiber.  Once  the  squeeze  force  on  the  i=10  fiber  begins  to 
increase,  the  normal  force  must  decrease  and  the  velocity  of  the 
i=10  fiber  decreases  from  the  free-fall  velocity.  Once  the  squeeze 
gap  reaches  a  value  of  0.24  <  at  t&3.2  seconds) ,  the  squeeze  force 
and  normal  force  acting  on  the  1=10  fiber  are  equal.  Since  the 
1=11,12  consolidation  fibers  are  static,  the  interaction  o+  the 
i=10  fiber  and  the  front  does  not  produce  any  additional  motion  of 
the  front.  Thus,  the  squeeze  gap  between  the  1=10  fiber  and  the 
front  continues  to  decrease,  until  the  i=10  fiber  reaches  the 
consolidation  front.  The  i  =10-12'  fibers  form  the  new  static 
consol i dat i on  front,  and  the  i =9  fiber  moves  toward  it.  This 


process  continues  until  all  12  or  n  fibers  are  stacked. 
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When  the  normal  gap  is  reduced,  the  -free-fall  velocity  of  the 
fibers  is  aecreased  according  to  equation  (71).  This  reduction 
again  arises  from  the  increased  resistance  of  the  normal  gap  to 
flow  through  the  more  narrow  channel.  The  consolidation  is  much 
slower,  and  the  sice  of  the  squeeze  gap  for  which  the  squeeze  force 
and  normal  force  are  equal  is  decreased.  From  equation  (70),  with 
£»o=0.5R,  the  value  of  the  "critical  squeeze  gap"  is  reduced  to  .11, 
and  the  free-fall  velocity  is  reduced  by  a  factor  of  3.2.  Halt  the 
dimensionless  distance  between  fiber  rows,  the  dimensionless 
velocities  and  the  dimensionless  forces  for  b  =0.5R  are  shown  in 
figure  7(a-c).  The  fiber  velocities  are  approaching  step 
functions,  indicating  that  the  squeeze  force  remains  small  for  a 
longer  period  of  time,  and  that  the  size  of  the  squeeze  gap  only 
contributes  to  the  force  when  the  gap  is  extremely  small.  As  the 
normal  gap  is  reduce  further,  the  normal  force  will  remain 
dominant,  due  to  the  large  resistance  to  the  normal  flow.  Only 
when  the  squeeze  gap  approaches  zero  will  the  squeeze  force  rapidly 
increase  and  slow  the  fiber  from  the  free-fall  velocity. 

5.  CONCLUSIONS  AND  DISCUSSION 

This  paper  treats  the  vertical  motion  of  long  fibers  through 
a  Newtonian  resin.  The  horizontal  distance,  260,  between  adjacent 
fiber  columns  remains  fixed,  while  the  vertical  distance,  21v  <t), 
between  the  adjacent  fiber  rows  i  and  ) + 1  decreases  with  time.  The 


fiber  motion  is  driven  by:  (1)  a  force  on  the  top  fiber  row  due  to 
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an  external  consol  i  dat  1  on  -force,  or  (2)  a  force  on  each  fiber  due 
to  its  tension  and  the  curvature  in  a  filament-wound  cylinder,  or 
(3)  any  combination  of  these  two  forces.  Lindt  C13  treats  the 
■first  case  with  only  a  force  on  the  top  fiber  row.  Lindt 
simplifies  the  two-dimensional  inertialess  fluid  motion  equations 
"using  a  modified  lubrication  approximation".  We  present  a 
rigorous  asymptotic  solution  for  c  =  ( £>0/R)  *  /  1  and 
£  =  ( h  /R )  * /'i'  <<i  ,  where  R  is  the  fiber  radius.  Unfortunately,  Lindt 
does  not  present  any  formulae  for  the  fiber  velocities,  pressures 
or  forces,  so  that  we  cannot  determine  the  difference  between  his 
lubrication  approximation  and  our  asymptotic  analysis. 

While  our  analysis  is  only  valid  for  e  < < 1  and  e . < < 1 ,  we 
present  results  for  e^l  at  t=0  and  en=l  in  order  to  compare  our 
results  in  figures  3  and  4  to  the  results  presented  by  Lindt  ill  in 
his  figures  6~lU,  for  the  same  values  of  h^(t-'O),  &Q,  R,  and  p. 
At  t=0,  Lindt  predicts  a  top  layer  velocity  in  excess  of  lmm/s, 
while  our  result  in  figure  3(b)  translates  to  a  dimensional  initial 
top  layer  velocity  which  is  two  orders  of  magnitude  smaller.  We 
suspect  that  there  is  an  error  in  the  scale?  for  Lindt  s  figure  8, 
as  a  top  layer  velocity  of  lmm/s  would  consolidate  the  top  two  rows 
in  0.02  seconds,  while  Lindt  indicates  that  this  consolidation 
takes  more  than  10  seconds.  A  more  fundamental  difference  concerns 
the  period  of  time  for  adjacent  fiber  rows  to  reach  equal 
velocities.  Lindt  predicts  that  the  top  two  rows  reach  an  equal 
velocity  at  roughly  7  seconds,  and  that  the  incremental  period  for 
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subsequent  rows  to  reach  equal  velocities  decreases  with  time.  The 
results  in  our  -figure  3(b)  indicate  that  the  top  two  -fiber  rows 

reach  equal  velocities  at  t=2.4  seconds,  and  that  the  period  tot- 

subsequent  row  consolidation  increases  with  time,  with  the  increase 
proportional  to  the  number  of  fiber  rows  consolidated.  Clearly 
there  are  some  significant  difference  between  the  lubrication 
solution  of  Lindt  LI]  and  the  present  rigorous  asymptotic  solution 
for  e  <<1  and  e 1  1 .  In  addition  to  yielding  the  correct  first 
order  solution  for  high  fiber  volume  fractions,  the  asymptotic 
analysis  defines  the  magnitudes  of  the  neglected  terms. 

The  results  for  a  force  on  all  of  the  fibers  show  that  the 

consol  1 dat i on  begins  at  the  bottom  fiber  row,  with  the  outer  fiber 

rows  moving  toward  the  compacted  fibers  in  the  interior.  We  have 

used  only  eleven  moving  fiber  rows.  However ,  even  this  small 

number  is  not  necessary,  since  the  development  of  distinct  wave 

fronts  would  permit  a  much  simplier  treatment.  The  wave  fronts  are 

more  abrupt  for  t  =(0.5)  ^ *  than  for  e  =■  i  .  As  c  and  e  decrease, 

n  n  n  i 

the  normal  flows  between  the  fiber  columns  becomes  progressively 
more  dominant.  Thus,  the  squeeze  flows  provide  negligible  or 
constant  forces  except  during  the  brief  periods  when  a  moving  fiber 
comes  very  close  to  a  static  fiber-. 

The  results  presented  here  and  by  Lindt  LID  provide  physical 
insights  into  the  local  characteristics  of  the  flows  during 
consolidation.  However,  for  a  Newtonian  resin,  the  macroscopic 
character i st l cs  of  consolidation  can  be  derived  from  D^rcy  s  law 
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■for  flow  through  a  porous  medium,  with  local  corrections  of  the 
permeability  as  the  local  fiber  volume  fraction  increases  L3J.  The 
more  important  value  of  the  present  model  is  its  potential  for 
treating  fiber-resin  motions  with  strongly  non-Newtonian  resins 
such  as  thermoplastics.  Darcy's  law  is  restricted  to  Newtonian 
fluids.  The  decoupling  of  the  normal  and  squeeze  flows  will 
provide  a  significant  simplification  for  a  non-Newtonian  resin. 
The  purpose  of  the  present  paper  is  to  provide  a  rigorous 
asymptotic  deri vatic  i  for  Lindt  's  lubrication  solution.  The 
accurate  definitions  of  relative  magnitudes  will  be  even  more 
important  for  a  non-Newtonian  fluid,  as  the  model  with  a  non- 
Newtom  an  fluid  will  involve  additional  dimensionless  parameters, 
time  scales,  etc.  In  a  future  paper,  we  will  present  extensions  of 
the  present  analysis  to  some  non- Newton i an  constitutive  equations 
for  thermoset  and  thermoplastic  resins. 
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FIG.  1  The  geometry  of  the  fiber-resin  system,  with  tne  normal 
and  the  squeeze  flow  regions. 

Fig.  2  Close-up  of  the  flow  regions,  (a)  the  squeeze  flow 
region  with  adjacent  fiber  row?  moving  with  velocities  V.  (top  row) 
and  LL  +  ^  (bottom  row.  The  rows  are  separated  by  a  distance  2hi  (t), 
and  the  fiber  surface  is  y,-,— f .  (b)  The  normal  flow  region 
between  adjacent  columns  separated  by  a  constant  gap  size  26.  The 
fiber  surface  is  at  x=f  <y),  and  Q  indicates  the  volume  flu;:  per 
unit  length  through  the  ith  normal  gap. 

Fig.  3  Results  for  the  Lindt's  problem  of  uniform  top  loading 
of  the  i=l  fiber  row  only.  The  iorcs  per  unit  lenth  is  T  =0.005 
N/m,  p=10kg/m-s,  fiber  radius  R=10  Jm,  and  h x ( t=0) = &Q=R.  (a>  A  plot 
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5.  Note  the  continually  decreasing  free  fall  velocity. 

Fiq.  4  Force  results  for  Lindt's  problem  o+  uniform  top  loading 
of  the  i  =  l  fiber  row  only,  (a)  The  forces  on  the  i=l  fiber 
separated  into  contributions  from  the  squeezing  flow  and  normal 
flow,  (b)  The  forces  on  the  l -2  to  4  fiber  rows  separated  into  the 
contributions  from  the  squeezing  and  normal  flows.  The  force  due 
to  the  normal  flow  is  greater  than  or  equal  to  zero,  wfu  le  the 
force  due  to  the  squeezing  flow  is  less  than  or  equal  to  zero. 
From  left  to  right  are  the  1-2,3, 4  force  contribution. 


Fig. 5  Results  for  Lindt's  problem  of  uniform  top  loading  of 

the  i=l  fiber  row  only,  but  with  a  reduced  normal  gap  size, 
&o=0.5R.  (a)  The  dimensionless  fiber  velocities  VL  for  i=l  to  3. 

Comparison  with  the  larger  normal  gap  results  show  that  the  free 
fall  velocity  for  the  consol i dat i on  front  is  reduced.  (b)  The 
forces  on  the  i=l  fiber  due  to  the  squeezing  and  normal  flows.  'cl 

the  forces  on  the  l  =2 , 3  fiber  rows  due  to  the  squeezing  and  nor  dial 

flows.  From  the  left  the  forces  on  the  i=2  fiber  row  are  shown, 

with  the  normal  force  on  the  i =3  fiber  increasing  from  zero,  the 
squeeze  force  decreasing  from  zero,  at  t=5  seconds. 

Fig.  6  Results  for  the  uniform  fiber  tension  on  all  l =n  fiber 

rows.  The  force  per  unit  length  is  T  =0.005N/m,  with  p= 1 Okg /m-s , 

and  h i < t=0 ) =6o=R= 1 0  um.  The  i=n  fiber  row  is  considered 

stationary.  (a)  Half  the  dimensionless  distance  between  adjacebt 

fiber  rows,  h  ^ ,  versus  time,  for  l =7  to  11.  Note  that 

consolidation  starts  at  the  i=ll  row  and  progresses  upward.  <b) 

*  _ 

Dimensionless  fiber  velocity,  V  ,  versus  time,  for  i=/  to  11  fiber 
rows.  Note  that  all  of  the  fibers  start  with  the  free  fall 
velocity  given  in  equation  (71)  with  T ^  replaced  by  T  .  (c)  the 

forces  on  the  fibers  separated  into  contributions  from  the 
squeezing  flow  and  normal  flow.  For  the  case  of  uniform  tension  on 
all  the  fibers  rows,  the  rows  which  have  consolidated  are 
stationary,  such  that  V1+1=0.  thus,  tne  normal  force  on  each  fiber 
row  decreases  from  one  to  zero,  while  the  squeeze  force  increases 
from  zero  to  one  during  the  time  each  fiber  is  moving  downward. 
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Fig.  7  Results  -for  uniform  tesion  on  all  i  =n  -fiber  rows,  but 
with  a  reduced  normal  gap  size  &  =0.5R.  (a)  Half  the  dimensionless 
distance  between  adjacent  fiber  rows,  h^ ,  versus  time  for  i = 1 0 , 11- 
This  case  conforms  more  closely  to  the  ideal  free  fall  behavior, 
(b)  The  dimensionless  fiber  velocities,  V*  ,  versus  time  for  i =9  to 
11.  The  initial  free  fall  velocity  is  reduced  due  to  the  decrease 
in  the  normal  gap  size.  (c)  The  forces  on  the  i = 1 1  fiber. 
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r'i 

lc 

eta  sub  "eye" 

1 

’  V  1  * 

*  * 

*1’  Vi 

lc 

and  y,  sub  one. 

sub  one  super  * 

X  -5 

X- 

-  y2» 

*  * 

*2’  V2 

1  c 

x  and  y,  sub  two. 

sub  two  super  * 

X  n 

y,  k 

*  * 
i  y 

lc 

x  and  y,  super  * 

3  partial  derivative 

Laplacian  operator 


V 


P  1  .  P* 

* 

P2>  P2 

pn »  Pn 
t,  t* 

At  ,  At  * 


u„ 

1 


gradient  operator 

ic:  p,  sub  one,  sub  one  super  * 

lc  p,  sub  two,  sub  two  super  * 

1 c  p ,  sub  1 c  n ,  sub  1 c  n  super  * 

ic  "tee",  super  * 

cap  delta  lc  "tee",  super  * 


1 


lc  "ewe"  sub  lc 


sub  sub  one 


4 


as 

above 

with  super  * 

uVi, 

* 

LVi 

lc 

"ewe” 

sub  1 c  y  sub 

sub  one 

as 

above 

with  super  * 

ux^’ 

* 

Uw 

1  c 

li  _  ,  _ _  ii 

ewt 

sub  lc  x  sub 

sub  two 

jC. 

aL 

as 

above 

with  super  * 

uy2’ 

* 

uy  2 

lc 

"ewe  " 

sub  1 c  y  sub 

sub  two 

as 

above 

with  super  # 

U  v  , 

u* 

lc 

"  ewe  " 

sub  lc  x  sub 

sub  1 c  n 

“n 

x  n 

as 

above 

with  super  * 

<-iv  , 

7  n 

* 

Uy 

7  n 

lc 

"  ewe  " 

sub  1 c  y  sub 

sub  1  c  n 

as 

above 

with  super  * 

U;. 

cap 

"ewe 

"  sub  1 c  x 

LL 

> 

cap 

"  vee 

"  ,  cap  "pee" 

£i'  £ 

n 

1  c 

epsilon  sub  lc  "eye",  sub  lc  i 

i  2  < ■<  2 

j  ,  f  2 

1  c 

-f  sub 

two  o-f  lc 

sub  two. 

1  c 

■f  sub 

two 

^n(y> 

lc 

+  sub 

len  o-f  lc  y, 

lc  -f  sub  1 

ei 

cap 

"cue 

"  sub  lc  "eye 

1  c 

"cue" 

sub  lc  "eye" 

Fext 

c  ap 

f  super  lc  "eye" 

sub  lc  "ext 

FJ 

c  ap 

F  super  lc  "eye" 

sub  lc  "jay 

<4 

1  c 

si  quia 

sub  lc  "jay" 

"kay" 

super  1 c 

"  eye  " 

n 

1  c 

n  (vector) 

nL, 

lc 

n  sub 

lc"kay" 

44 


W 


r 


cap  "tee"  sub  lc  "eye",  super  * 
1c  theta 

cap  delta  lc  "pee"  sub  lc  n 
cap  "tee"  sub  lc  "see" 
lc  "tee"  sub  lc  "see" 
scri pt  1 c  "el " 
cap  "double  ewe"  sub  lc  r 


cap  F  super  lc  "eye"  sub  lc 
sub  sub  lc  "sq" 

cap  F  super  lc  "eye"  sub  lc  y 
sub  sub  lc  "sq" 

cap  F  super  lc  "eye"  sub  lc  y 
sub  sub  lc  "sq,bot" 

cap  F  super  lc  "eye"  sub  lc  y 
sub  sub  It:  "sq,top" 

cap  F  super  lc.  "eye"  sub  lc 
sub  sub  1 c  n 


cap  F  super  lc  "eye"  sub  lc  y 
sub  sub  lc  n 


TIME  (SECONDS) 


Fl§  bCOk)  T 


TIME  (SECONDS) 


VOG&Kof 


TIME  (SECONDS) 


HQ 


